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The search for neutral resonances at the energy frontier has a long and illustrious history, resulting 
. in multiple discoveries. The canonical search scans the reconstructed invariant mass distribution of 

' identified fermion pairs. Two recent analyses from the CDF experiment at the Fermilab Tevatron 

have applied novel methods to resonance searches. One analysis uses simulated templates to fit 
(3\ , the inverse mass distribution of muon pairs, a quantity with approximately constant resolution for 

momenta measured with a tracking detector. The other analysis measures the angular distribution of 
electron pairs as a function of dielectron mass, gaining sensitivity over a probe of the mass spectrum 
alone. After reviewing several models that predict new neutral resonances, we discuss these CDF 
analyses and potential future applications. 



I. INTRODUCTION 



Searches for neutral resonances have historically brought major breakthroughs by either confirming important 
predictions or discovering unexpected particles. The 1974 discovery of the J/ip meson [l| as a cc bound state confirmed 
the GIM mechanism [2j for preventing flavor-changing neutral currents, and the discovery of the Z boson [3] confirmed 
the gauge unification of the electromagnetic and weak forces Meanwhile, the discovery of the upsilon [f| was 
completely unexpected, and increased the number of known fermion generations to three. 

Turning to the future, there are reasons to expect the next important particle physics discovery will be a neutral 
resonance. In addition to the well-motivated Higgs boson Q of the standard model (SM), there are many new 
| resonances predicted by proposed extensions to the standard model. These extended theories can address unexplained 
features of the SM, such as: the lack of gauge unification and the hierarchy between the electroweak and Planck scales 
(through supersymmetry Q or the presence of extra dimensions Q); and parity violation and light neutrino masses 
• i-h , (through an additional SU(2)r gauge symmetry [jj], which has weak couplings to right-handed fermions). 

The most sensitive direct searches for neutral resonances at high mass come from Tevatron pp collision data. Future 
searches in pp collisions from the Large Hadron Collider (LHC) [l(| will increase the probed mass range. As larger 
datasets with higher energies are studied, enhancements to the search strategy can improve sensitivity and facilitate 
the analysis. Such enhancements have been developed and applied to searches for new resonances in the CDF dimuon 
[ill and dielectron [L21 data. 



II. MODELS CONTAINING NEUTRAL RESONANCES 



A neutral resonance decaying to fermion pairs can have intrinsic spin equal to 0, 1, or 2. No fundamental scalar 
particle has yet been observed, though the SM requires one in the form of a Higgs boson. Beyond the SM, there 
could be multiple Higgs bosons with varying properties fl3| . In supersymrnetric models, there are spin-0 partners to 
fermions that could be produced as resonances in pp or pp collisions [14j |. Any model with an additional U(l) gauge 
group will have a new spin-1 gauge boson, generically referred to as a Z' boson [15|, [l6|. Models of extra dimensions 
at the electroweak scale predict spin- 2 graviton resonances [l7l. [l8l|. 
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A. Sneutrino Production in Hadron Collisions 



To remove the fine-tuning of the Higgs boson mass, the scale of supersymmetry should be of the same order as 
the electroweak scale, making the discovery of supersymmetry likely at the Tevatron or LHC (should it exist). The 
supersymmetric partner to the neutrino has no electromagnetic charge and is thus a candidate for production as a 
neutral resonance. 

Resonant sneutrino production would violate i?-parity, a multiplicative quantum number that is +1 for matter and 
-1 for supersymmetric matter. The violation of i?-parity implies that the lightest sparticle is not stable, potentially 
removing it as a candidate for dark matter. However, if the lightest sparticle has sufficiently small couplings to give 
it a lifetime on the order of the age of the universe, it can still be a dark-matter candidate 

Proton decay limits require at least one set of i?-parity- violating terms to be vanishingly small [201 ] . This can 
be accomplished by imposing a "baryon parity" that conserves baryon number and suppresses proton decay more 
than i?-parity conservation [2l|. With baryon parity there are two sets of i?-parity- violating Yukawa terms in the 
Lagrangian, both of which are consistent with proton decay limits and relevant for sneutrino production and decay 
at a hadron collider: 



£if. - /.,/. ,< \ + X ljk L t Q 3 d c k , (1) 

where Q (d) and L (e) are SU(2)l doublet (singlet) superfields. As shown in Fig.[l] the first term governs sneutrino 
decay to leptons and the second term governs sneutrino production in hadron collisions. Because of the d k superfield 
in the production term, only-down type quark interactions produce sneutrinos, with \\\\ the most relevant coupling 
at the Tevatron and LHC (due to parton distributions in the proton). 




FIG. 1: The Feynman diagram for resonant sneutrino production at a hadron collider. 



The total width of a new resonance is an important parameter in a search, in particular relative to the detector 
mass resolution. The partial width of a given sneutrino decay is [22I [23| : 

m^fjfk) = ^X 2 m Di , (2) 

where Cj k is a color factor, and A (— \j k or Ay fe ) is the coupling to the final-state fjfk- The width is fairly narrow; 
for example, if only X' al and one Xijk are large, then for respective values of 1/2 and 1, the width is 3.5% of the mass 
(this can be compared to the Z boson, whose width is 2.8% of its mass 24]). 

A range of indirect limits exist on \' ixl 25], and generally depend on other supersymmetric parameters. A typical 
set of limits comes from the ratio of r(7r — > ev)/Y{Ti — > fiv) or T(r — > niy)/T(Tr — > [iv)\ A^ xl < dim^/TeV, where 
a\ = 0.26, 02 = 0.59, and 03 = 1.2. For reasonable values of mg (0.2 — 1 TeV), these limits allow significant production 
rates with a relatively narrow decay width. 



B. Z' Vector Bosons 



Many models predict a new electroweak-scale U(l) gauge symmetry [26|, [27], HH, Hil , which would have an associated 
Z' gauge boson. A useful test model is a superstring- inspired grand unified theory with EgxE 8 gauge structure (30j . In 
this model, the Eg group is a hidden sector that breaks supersymmetry, and the Eg group is broken to Eg x SU(3) by the 
compactification of extra dimensions [oil ]. Each generation of matter particles fits in a fundamental 27 representation 
of E 6 ; thus, before E 6 is broken, each SM generation is just a single field distinguished by its E 6 charge. 
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The range of options for breaking Eg to the SM gauge structure allows a variety of phenomena [16| . A symmetry- 
breaking proceeding through SO(10)xU(l) to SU(4)cxSU(2) l xSU(2) r xU(1)b-l restores parity conservation and 
provides for the seesaw mechanism for small neutrino masses. Alternatively, the breaking can proceed through 
SO(10)xU(l)^, to SU(5)xU(l) x xU(l)^,, producing two new U(l) gauge groups. At a high mass scale, the SU(5) can 
be broken to the SM gauge groups and one of the extra U(l) gauge groups can be broken, potentially leaving one 
non-SM U(l)' at the electroweak scale. Taking this U(l)' to be a linear combination of U(l),/, and U(l) x , 

U(l)' =U(1)^ cosS + U(1) X sin0, (3) 

a generic U(l)' can be expressed in terms of 8 [32j. 

Scanning the 8 parameter space gives models with distinct phenomena. The secluded U(l)' (8 = it — tan -1 y / 27/5) 
is mediated by a Z' sec boson whose mass results from the vacuum expectation value of a scalar held with no 
SM charge [Hj]. For 8 = — tan -1 yT/15, the right-handed neutrino has no charge in the extended gauge group 
(SU(3) c xSU(2) L xU(l) Y xU(l) N ), and is thus sterile [H]. The breaking of E 6 directly to the SM groups plus U(l), ; 
corresponds to 8 = tan -1 \/3/5. If the breaking proceeds through an extra SU(2)i group (instead of SU(2)r), then 
the W and Z 1 bosons of the new group have zero electromagnetic charge. 

More general classes of models have also been considered, with the constraint of anomaly cancellation to produce 
a consistent theory [35]. Under the assumption of the SM Higgs mechanism for generating fermion masses, a general 
class of models has U(l)' charge B — xL, where B (L) is baryon (lepton) number (and the right-handed neutrino 
charge is fixed to -1). Allowing for non-SM mass generation but considering only SM particles for anomaly cancellation 
gives U(l)' charges 1/3, x/3, (2 — x)/3, -1, —(2 + x)/3, and (—4 + x)/3 for the states qx, ur, dp, II, e_R, and vr, 
respectively. This model is referred to as q + xu and includes the case of B — L symmetry for x = 1 . Two additional 
model classes arise when two non-SM fermions are added to the theory. One is referred to as d — xu and has charges 
of 0, 1/3, and —x/3 for q^, da, and ur, respectively. The other, 10 + x5, has fermions in the 10 representation of the 
SU(5) grand unified group with the U{1)' charge 1/3, and fermions in the 5 representation with charge x/3. 

In general, couplings of the new Z' boson to SM particles are smaller than those of the Z boson in the SM. However, 
the new Z' boson could decay into the non-SM particles that are part of the 27 representation of Eg. If decay to all 
of these particles is possible, the Z' boson width could be «5% of its mass 36]. Even in this extreme case, the Z' 
boson would appear as a narrow resonance. 

C. Graviton Resonances 

It has been suggested that the apparent difference between the scales of g ravity and electroweak symmetry-breaking 
is due to the presence of at least one unobserved spatial dimension [1, flTl ]. The spread of the gravitational field into 
the extra dimension(s) weakens the strength of gravity in the observed dimensions. Closing the gap between the 
electroweak and Planck scales requires either the number or the size of the extra dimensions to be large, if they are 
flat. 

Recently, Randall and Sundrum have proposed a model that removes the scale hierarchy using one small extra 
dimension [17] . This can be accomplished with a warped dimension separating the SM brane from the gravity brane, 
resulting in a metric of the form: 

ds 2 = e-^^dx^dx" - r 2 d(f> 2 , (4) 

where r is the compactification radius, and k 2 and <ft = [ — ^] are the spacetime curvature and coordinate in the 

extra dimension, respectively. The curvature is of the order of M Pl , where Mpi — Gjy 1 ^ 2 ~ 10 19 GeV is the Planck 
scale on the 4-dimensional spacetime and Gat is the Newtonian gravitational constant. For a string theory with 0(1) 
couplings, k/Mpi ~ 0.01 QJ]. 

Due to the spacetime warping, distances are exponentially larger on the gravity brane, resulting in a large gravita- 
tional field flux on this brane. The gravitational force for an observer on the SM brane appears as: 

F~m im2 /(M 2 EW e 2kr7T R 2 ), (5) 

where R is the distance between masses mi and ni2 in the three large spatial dimensions. Thus, kr ~ 12 reproduces the 
observed weakness of gravity and there are no large hierarchies in the model. In terms of the gravitational quantum, 
the wave function of the massless graviton state is localized on the gravity brane and exponentially suppressed on the 
SM brane. 
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Graviton excitations are localized on the SM brane and are thus expected to have mases m n at the electroweak 
scale, 

m n = kx n e- krv , (6) 

where x n are 0(1) roots of a Bessel function and ke~ kr7T is of the order of the electroweak scale. The resonance width 
is proportional to (k/Mpi) 2 , and is less than a few percent for k/Mpi < 0.1. 



III. COLLIDER SEARCHES FOR NEUTRAL RESONANCES 



The most stringent direct limits on new neutral resonances come from searches at the Tevatron. Run II searches 
at the CDF and DO experiments have probed resonance decays to 
"MElHIiil, light quarks [HHl, top quarks 1 



mirs of electrons 




taus 
and Z 



47(, gluons [4J, |46j, photons |4j 




50] bosons. The most sensitive searches use long-lived final-state particles, while other searches 
cover parameter space where couplings to long-lived particles are suppressed. 

Searches for a resonance decaying to a pair of stable particles typically probe the reconstructed invariant mass 
distribution for evidence of a narrow peak, with the peak width determined primarily by detector resolution. Because 
detector resolutions increase with increasing mass, the expected peak width also increases. This complication usually 
results in search windows that change as a function of mass, with the window causing some loss of acceptance. A 
recent CDF search [ill ] uses templates to fit the full dimuon inverse invariant mass spectrum for new resonances, 
avoiding the acceptance loss from a search window. In addition, the inverse mass distribution has approximately 
constant resolution, simplifying the search. 

A complement to the invariant mass distribution is the angular distribution of the final-state particles, which can 
be used to separate a signal from the SM background and to determine the spin of the new resonance. CDF has 
performed a search in the dielectron final state using the cos (9* distribution [12], where 8* is the angle between the 
electron and the incoming quark in the boson rest frame [5l| . 



A. CDF Dimuon Search 



Currently, the CDF dimuon analysis of 2.3 fb _1 of pp collision data is the most sensitive search for neutral resonances 
over most of the probed parameter space. The search uses a parametric simulation to model the detector response 
and resolution for muons, calibrated using known resonances. After normalizing the SM inverse mass spectrum to 
the Z boson peak, the data are fit as functions of the number of new-neutral-resonance events above background and 
the resonance pole mass. No statistically significant excess above the background is observed, and limits are set for 
the various test models. 



1. Detector Alignment and Calibration 

In the CDF search, muon momenta are measured with the central outer tracker (COT) [52 ]. a wire drift chamber 
embedded in a 1.4 T magnetic field covering \rj\ < 1 and radii 43 cm to 133 cm 53:] ■ The reconstructed tracks 
are constrained to originate from the time-averaged transverse beam collision coordinate, significantly improving 
momentum measurement resolution. Calorimeters and a muon detector system at large radii from the beam line are 
used for muon identification and event triggering. 

To minimize bias and optimize the detector resolution (and thus the statistical significance of a narrow resonance), a 
detailed alignment of the COT is performed [5~ij . The alignment uses cosmic-ray muons reconstructed as a single track 
through both sides of the nominal collision point [55]. With 96 radial wire layers in the COT, any given two-sided 
track has up to 192 measurement points. The 96 layers are divided into 8 superlayers of 12 wires each, with each 
superlayer containing enough 2 cm wide cells to cover the azimuth. The first stage of the alignment allows a rotation 
and a shift of each cell, such that the mean residual of hits in any given cell is statistically consistent with zero (with 
a precision of a few microns). 

After the individual cell alignment, a global correction to the wire shape between endplates is derived as a function of 
(j) and radius. The shape is determined by the gravitational sag from the weight of the wire, and by the electrostatic 
deflection from the local electric field. The nominal correction to the wire shape due to these effects is further 
modified with an empirical correction function derived from measured biases between the two separately fit sides of 
the cosmic-ray tracks. 
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A final correction to the track curvature is applied after the track reconstruction. The correction is derived from 
the difference in the ratio of calorimeter energy to track momentum for positrons and electrons, as functions of <fi and 
cot 

The momentum scale is calibrated by tuning the measured J/tp, T > an d Z boson masses to their precisely known 
values. Individual hit resolutions of 150 microns are determined from the observed width of the T — * peak, 
consistent with hit residuals of muons from Z boson decays. The transverse beam profile is modelled as a gaussian 
with a size set by the observed width of the Z — > mi mass peak measured with beam-constrained tracks. 

2. Inverse Mass Scan 

Muon momenta transverse to the beam line are determined from a measurement of the reconstructed track curvature 
c. The Lorentz force qv x B causes a helical trajectory of the muon, resulting in a transverse momentum: 

p T = eBR, (7) 

where R is the radius of curvature and B is the magnitude of the magnetic field. 

The resolution on the track curvature (c = i? _1 /2 at CDF) can be derived from the individual hit resolutions. 
Taking as an example a muon produced with no impact parameter at (j> = 0, its position at a given radius r is 
(rcos<fi,rsm(f>). Using the equation of the track circle, x 2 + (y — R) 2 = R 2 , the curvature is c = smcf>/r and its 
resolution due to a measured hit is 5c = cos 4>5<j)/r. The partial derivative of the curvature with respect to the hit 
resolution 5D = rd(f> is then 

Sc/SD = (1 - c 2 r 2 f/ 2 /r 2 . (8) 

For muons with small curvature (large momentum) , the resolution is effectively independent of curvature and improves 
with the square of the detector radius. 

Defining the reconstructed muon energies to be E\ and -E2, and their opening angle to be 0, the measured mass 
can be expressed as: 

m=[2E 1 E 2 (l~cosQ)] 1/2 . (9) 

The dominant contribution to the mass resolution comes from the momentum measurement, since the angular resolu- 
tion is negligible by comparison. A high-mass resonance is predominantly produced with a relatively small transverse 
boost, resulting in muons with similar transverse momenta. Then, m a pt, or 1/m oc c. Thus, a new narrow 
resonance would have an approximately constant width in the reconstructed 1/m distribution of central muons. 

The search in a constant-width variable simplifies the analysis. A distribution with uniform binning can be visually 
scanned for resonances. Taking bin widths sufficiently narrow with respect to the resolution allows a template fit for 
a resonance centered on each bin, with the step size an equal fraction of the peak width throughout the distribution. 
This procedure optimizes the scan of the 1 /m distribution. 

At CDF, the inverse mass resolution is 17% TeV -1 , with an additional contribution from multiple scattering at 
low mass. At 100 GeV, the total resolution is about 30% smaller than the intrinsic width of the Z boson. Thus, the 
width of the new resonance could noticeably broaden the peak at low mass. However, the relative resolution increases 
linearly with increasing mass, while the relative instrinsic width remains constant, so for most resonances the detector 
resolution will dominate above a few hundred GeV. 

In the CDF analysis, the search region is 35 bins of m^ 1 < 10 TeV -1 , resulting in an expected peak width of about 3 
bins due to detector resolution. The 70-100 GeV mass range is used to normalize the expected background, effectively 
removing systematic uncertainties due to luminosity, background cross section, and trigger and muon identification 
efficiencies. For each probed inverse mass bin, a template of combined signal and background is compared to the data 
to determine the number of signal events that maximizes the log Poisson likelihood. 

The template fit adds acceptance outside of the usual mass window, particularly at masses near the kinematic 
threshold. While the template neglects interference between the new resonance and Z and 7 bosons, any interference 
has a small effect on the search because the resonance is narrow and has a small cross section [35j ]. 

3. Backgrounds 

Dimuon production at a hadron collider occurs predominantly through the Drell-Yan process of Z/7* production. 
This process has been calculated at next-to-next-to-leading order (NNLO) in a s [3|| and next-to-leading order (NLO) 
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in ctEW [56| at the Tevatron, and is thus well understood theoretically. Since the CDF analysis normalizes the 
background to the Z boson peak, the shape of the inverse mass distribution is the important theoretical input. Over 
the probed mass region, the ratio of NNLO to leading order (LO) predictions have an w 10% variation. The difference 
between NLO and NNLO predictions is taken as a systematic uncertainty and increases from zero at 91 GeV (the 
normalization mass) to 9% for a mass of 1 TeV. A more important systematic uncertainty at 1 TeV comes from 
uncertainties on the parton distribution functions (PDFs). At this mass, the valence (anti-)quarks must have a large 
fraction of the (anti-)proton's momentum. The uncertainty obtained from a comparison of CTEQ [57] PDFs fit using 
LO and NLO inputs increases with mass to approximately 16% at 1 TeV. No NLO chew correction is applied, and a 
3% uncertainty at 1 TeV is incorporated to cover its neglect. 

A background relevant at high mass arises from hadrons decaying to muons in the COT, where the track is mis- 
reconstructed as a nearly straight line. The misreconstruction occurs via a track kink at the decay vertex or other 
incorrectly assigned hits in the inner superlayers. This background is reduced by requiring each muon's COT hit pat- 
tern and track fit y 2 to be consistent with a well-measured track resulting from the collision vertex. The small residual 
background is estimated using dimuons with the same charge, which are assumed to arise from misidcntification or 
other non-prompt sources. 

W boson pairs produced either directly or through top quark decays contribute to the dimuon sample when both 
W bosons decay to [iv. These backgrounds have been calculated to NLO in a s [H, [5!| and are only relevant at high 
mass. 

A final potential background arises from cosmic-ray muons passing through the detector. These muons can have 
high energies and thus contribute to the background at high mass. At CDF the cosmic-ray background is effectively 
eliminated by the two-sided track fit used to identify cosmic-ray muons in the COT [55| . 

4- Signal Cross Sections and Acceptance 

Cross sections for sneutrino production have been calculated at NLO in a s [6(3] • If there is CP conservation in the 
sneutrino sector, the anti-sneutrino and sneutrino cross sections will be the same; this assumption is made in the CDF 
analysis. Z'-boson and Randall- Sundrum-graviton cross sections are determined from pythia I6l| with an NNLO a s 
correction factor applied [35| , and the appropriate couplings are used for the various Eg models |62l] . 

In the CDF analysis, acceptance is calculated as a function of spin and inverse mass. For 1 TeV resonances, which 
are produced with relatively little longitudinal momenta, the acceptance for observing two central muons is ~ 40%. 
At 100 GeV, the larger average longitudinal boost reduces the acceptance for the CDF central muon and tracking 
chambers, resulting in an acceptance of ~ 15%. A relative uncertainty of 3% is estimated from comparisons of the 
parametric and GEANT-based detector simulations. A study of Z — ► /.ifi events shows the momentum dependence of 
the muon identification efficiency to be well modelled by the simulation. 

5. Results and Limits 

The observed dimuon invariant mass spectrum (Fig. ^ shows no evidence for a new neutral resonance. This can be 
clearly seen in the error- weighted difference between observation and expectation (Fig. [3]) . A narrow resonance would 
appear as a significant excess in one bin, with additional excesses in neighboring bins. The most significant observed 
excess occurs at the lowest probed mass in the search, 103 GeV. An ensemble of simulated experiments gives a 6.6% 
probability of finding a more significant excess anywhere in the search region from a background fluctuation. 

The results are translated into cross-section limits for new resonance production (Figs. 0][6]), and then into mass 
limits for specific models (Table JJ. The mass limits are the highest of any search, except for models with weak 
couplings. For such couplings, the probed mass range is lower and the recent CDF dielectron search [371 ] has better 
sensitivity, since the CDF calorimeter has broader coverage in rj. The lower the resonance mass, the larger the average 
boost in the beam direction, and the larger the average rapidities of the leptons. The dielectron analysis thus has 
better acceptance and sensitivity at these masses. 

B. CDF Dielectron Search 

The CDF dielectron search in 0.45 fb _1 of pp collision data is the only hadron-collider search for Z' bosons to 
use the dilepton angular information. The search parametrizes the detector response in the (m ee ,cos#*) plane and 
distinguishes the SM and Z'-boson hypotheses using the observed distribution in this plane. The data are consistent 
with the SM so limits are set on Z' bosons in a generalized model parameter space. 
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FIG. 2: The inverse invariant dimuon mass for data and expected background. The Z boson peak is prominent at » 11 TeV , 
and a new resonance would appear as a similar (narrower) peak in the sloping region on the left. 
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FIG. 3: The difference between data and expected background, divided by the statistical uncertainty. The data are consistent 
with the background expectation; the most significant excess appears in the rightmost bin, corresponding to a mass of 103 
GeV. 



1. Two-dimensional (m ee , cos 0*) Scan 



The scattering amplitude for the // — > e e + process is [6J 



Mi = Mfif -» eje+) 



+ ^ - 



-Qe 2 



§—m/%-\-imzT z 



Cf(f)C?(e) + 



Cf(f)Cf(e), 



Z' i 



s — m , -\-im < 



(10) 



where i and j are the fermion helicities (L, R), Q is the electromagnetic charge of fermion /, C i j (/) are the fermion 
couplings to the Z and Z 1 bosons, mz (£z) an d m z> (Pz 1 ) are the respective Z and Z' boson masses (widths), and 
s is the squared center-of-mass energy of the collision. Using this amplitude, the differential angular cross section is 



da 



1 



*\2l 



[(I All I + \Arb\ z ){1 +cos0*y + (\A LR \ 2 + \A RL \ 2 )(1 - cos9*) 



dcos9* 128?rs 

A Z' boson with non-zero couplings to quarks and electrons alters the SM cos 9* distribution. 



(11) 
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FIG. 4: The cross section limits for a new spin-0 resonance, and the theoretical predictions for sneutrino production for various 
values of the coupling squared (A^j = A 2 ) times the branching ratio to dimuons. 
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FIG. 5: The cross-section limits for a new spin-1 resonance, and the theoretical predictions for Z 1 bosons with the same 
couplings to fermions as the Z boson (Z' SM ), and in various Es-inspired models. 
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FIG. 6: The cross-section limits for a new spin-2 resonance, and the theoretical predictions for R-S gravitons for various values 
of k/Mpi. 
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TABLE I: 95% C.L. lower mass limits, in GeV, for sneutrinos, Z' bosons, and gravitons with various model parameters. The 
Z' SM boson has the same couplings to fermions as the Z boson. 



V 

(An) ■ 


V 

mass limit 


Zj 

model 


Z 

mass limit 


Kb graviton 


graviton 
mass limit 


0.0001 


397 


Z 1 , 
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0.01 
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0.0002 
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7' 


821 


0.015 


409 


0.0005 
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z N 
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0.025 


493 


0.001 


662 




878 


0.035 


651 


0.002 
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892 


0.05 


746 


0.005 


810 


Z v 


904 


0.07 


824 


0.01 


866 


Z'sm 


1030 


0.1 


921 



The CDF search probes the (m ee ,cos#*) plane in (10 GeV, 0.25) bins using a look-up table from the full detector 
simulation to determine the acceptance in each bin. The Z' boson signal hypothesis, including interference with the 
Drcll-Yan process, is compared to SM Z/j* boson production through a likelihood ratio for the two hypotheses in 
the search region of m ee > 200 GeV. The processes are modelled with pythia [6l[ and a mass-dependent NNLO 
correction factor [35j ] - 

In the search region CDF estimates the following SM backgrounds in 0.45 fb _1 of data: 80 Drell-Yan events; 28 
events with a jet misreconstructed as an electron (dijet and W + jet); and 7 diboson (WW and WZ) events. The 
misreconstructed-jet events are estimated by applying a jet-to-electron misreconstruction rate to all jets in events 
with one reconstructed electron and at least one jet. Diboson events are estimated with PYTHIA and their theoretical 
cross sections [58| . 



2. Results 



The 120 observed events in the search region are consistent with the HS^g expected background events in the 
two-dimensional (m ee ,cos0*) plane. The projections along the m ee and cos 9* axes are shown in Figs. [7] and [8] 




FIG. 7: The invariant dielectron mass for data and expected background. The inset shows the data events in the search region 
m ee > 200 GeV. 

From the comparison of SM and Z' hypotheses, limits are set on the masses and couplings of Z' bosons. The mass 
limits range from 675 GeV for Z' sec to 860 GeV for Z' SM . For the generalized models B ~ xL, q + xu, d — xu, and 
10 + x5, limits are set in the two-dimensional plane of (x,Mz'/gz>) for several values of gz 1 (Fig. [9]). The limits 
extend to smaller values of \x\ and gz' than those from LEP II. 
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IV. FUTURE SEARCHES 



The full Run II data sets are expected to be a factor of 4 larger than that used in the CDF dimuon analysis. 
The highest current mass limits are at the edge of a steeply falling parton-luminosity curve and will not significantly 
increase. However, there is ample opportunity to see hints of a weakly-coupled new resonance at a mass below the 
kinematic threshold using the full CDF and DO data sets. With the higher energy collisions soon expected from the 
LHC, the highest probed mass could increase by a factor of 3 or more [k}. Thus, neutral resonance searches will 
continue to provide significant discovery potential into the future. 

The technique applied to the CDF dimuon search can be extended to searches in other final states. For example, at 
high mass the calorimetric measurement of electrons and photons has a fractional energy resolution that is constant 
in energy because the calorimeter sampling term becomes negligible in comparison. Since 5\nm — 5m/ m, a new 
resonance decaying to electron pairs will have a constant peak width in the lnm distribution. This is not affected 
by the fractional intrinsic width, which is also constant in mass. For resonances decaying to quarks and gluons, the 
optimal distribution depends on the calorimeter. At low energy the resolution is proportional to V~E, so the fractional 
resolution improves with increasing energy. In this case resonances will have constant width in the \~E distribution. 
At sufficiently high mass, however, the intrinsic width and constant fractional resolution term will become dominant, 
in which case resonances will have constant width in lnm. 

The angular distributions of the final-state particles included in the Z'-boson search at CDF have also been used 
by DO to search for graviton production [64l loq l . Optimally, an unbinned likelihood can be performed for each spin 
hypothesis, using the full matrix-element and resolution information on an event-by-event basis. This technique has 
been applied to Higgs-boson and top-quark searches [49|, [66[ and measurements [13] at the Tevatron. For a high-mass 
resonance search, a modest gain in sensitivity is expected beyond a two-dimensional (mu, cos 9*) fit. 
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